VORTICES AND MAGNETIZATION IN KAC'S MODEL 
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Abstract. We consider a 2-dimensional planar rotator on a large, but finite lattice with 
' a ferromagnetic Kac potential J-yii) = 7^J(7z), J with compact support. The system is 
subject to boundary conditions with vorticity. Using a gradient-flow dynamics, we compute 
minimizers of the free energy functional at low temperature, i.e. in the regime of phase 



Vector spin models with an internal continuous symmetry group, such classical 0~^{q) 



■ transition. We have the numerical evidence of a vortex structure for minimizers, which 
> ■ 

CN ' present many common features with those of the Ginzburg-Landau functional. We extend 

^ I the results to spins valued in 5"^ and compare with the celebrated Belavin & Polyakov model. 

O ' 0. Introduction. 

o 

p ^. models (XY or "planar rotator" for q = 2, and Heisenberg model for q = 3,) play an important 
r-| ■ role in Statistical Physics. In one or two dimensions, and for all inverse temperature /3, if 
' the range of the translation invariant interaction is flnite, then a theorem of Dobrushin & 
Shlosman shows there is no breaking of the internal symmetry (that is, Gibbs states are 
invariant under 0~^{q)) and furthermore, by a theorem of Bricmont, Fontaine & Landau, 
uniqueness of the Gibbs state holds (see e.g. [Si,Chap.III]). 

Despite of this, a particular form for phase transition exists, which can be characterized 
by the change of behavior in the correlation functions. In the low temperature phase they 
have power law decay, showing that the system is in a long range order state (exhibiting in 
particular the so-called "spin waves", ) but they decay exponentially fast at high tempera- 
tures, breaking the long range order, even though thermodynamic quantities remain smooth 
across the transition. For the XY system, these transitions were described by Kosterlitz & 
Thouless in term of topological excitations called vortices : while these vortices are organized 
into dipoles at low temperature, a disordered state emerges at the transition. But the obser- 
vation of the spatial distribution of defects shows that it is not uniform ; rather, defects tend 
to cluster at temperatures slightly larger than the transition temperature, and there are still 
large ordered domains where the spins are almost parallel (see e.g. [LeVeRu], [BuPi],[MiZh], 
and references therein. ) 
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Here we consider a Kac version of the classical XY or Heisenberg model on a "large" 
lattice A C Z^. It was studied in particular by Butta & Picco [BuPi]. The hamiltonian 
(except for the interaction with the boundary) is of the form 



where 7 is a small coupling constant and J denotes a cutoff function. Kac potentials for fixed 
7 have finite interaction, but as we take an appropriate limit 7^0, they can be considered, 
to this respect, as long range. Thus, they share some features with the mean field model, 
though exhibiting better mechanisms of phase transitions, which depend in particular on the 
dimension, as for the short range case. For the mean field model with 0^{q) symmetry, 
= 2, 3, we know that there is no phase transition for inverse temperature P < q (Gibbs 
measure is supported at the absolute minimum of the free energy functional, ) while there is 
a phase transition for P > q, with internal symmetry group 0~^{q). 

When the model possesses internal symmetry and common features with the mean field, 
it is hard to expect vortices at low temperature, unless the symmetry is somehow broken, 
for instance if the system is subject to boundary conditions. This situation is met in other 
domains of Condensed Matter Physics, as in supraconductivity, where vorticity is created by 
an exterior magnetic fiux, or for superfluids. In that case, phase transitions of matter are well 
described by critical points of free energy (Ginzburg-Landau) functionals ([BeBrHe], [OvSi], 
etc. . . ) 

One of the main process consists in averaging the spins cta over some mesoscopic boxes, 
so to define the magnetization m = m\* on another "coarser" or "mesoscopic" lattice A*. 
The free energy (or excess free energy) functional Ffj^^{m) at inverse temperature f3 in case of 
Kac models with internal symmetry, can be simply derived from a suitable renormalization 
of making use of the entropy I{m) for the mean field that corresponds to Van der Waals 
free energy fp{m) = — i|mp + ^I{m) (see Sect.l). 

To understand the significance of Fp^~f{mA*), one should think also of the formal "sta- 
tionary phase" argument, as A ^ 00, which suggests that an important role in the averaging 
with respect to Gibbs measure, is played by configurations close to those which produce the 
local critical points of -Fg,7. This occurs in computing correlations functions (see e.g. [Z]. ) 
These critical points consist in ground states, or metastable states. 

They will be determined as the attractors of a certain dynamics, similar to this given by 
the "heat operator" , but known in that context as the gradient-fiow dynamics [DeMOrPrTr] , 
[DeM], [Pr] .... Thus, we expect convergence of this dynamics toward a Gibbsian equilibrium, 
though this will not be rigorously established here. 

Let us present our main results. 
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In Sect.l, we describe in detail Kac's hamiltonian on the lattice, and recall briefly the 
renormalization scheme, that makes of the free energy functional a fairly good approximation 
for the density of Gibbs measure, i.e. A*/3,7,a ~ exp[— /?7~^.F(m|m'^)]. Here J^{m\m^) denotes 
the free energy functional subject to boundary conditions on A**^. 

In Sect. 2, we present a simple, combinatorial averaging process, relating Kac's hamilto- 
nian with the free energy functional. While exhibiting the main idea of renormalization, it is 
more suitable for effective computations on the lattice. 

In Sect. 3, we study Euler-Lagrange equations for the free energy functional, and intro- 
duce the corresponding gradient-flow dynamics m{x, t) (see Eqn. (3.4). ) Using that J^{m\m'') 
is a Lyapunov function, we show that m{x, t) converges towards a critical point of T{m\m'^), 
generically, a local minimum. Unless /? < 2, in which case m = is the unique minimizer of 
T{m\m''), as expected from the considerations above on the mean field, in general there can- 
not be uniqueness of the limiting orbits, at least for a finite lattice. Instead, local minimizers 
might depend on initial conditions m{0,x) inside A*. 

Local minimizers however, have the property that their modulus be bounded by m/3, if 
this is true of the initial condition, and as expected from general results relative to the Gibbs 
states [BuPi], \m\ has to be close to mp on large regions of A*. Actually Proposition 3.4 
indicates that if no vorticity is induced by the boundary, nor by the initial condition, then 
all magnetizations of the limiting configuration should point out in the same direction and 
have length about mp. 

In Sect. 4 we make numerical simulations, introducing a boundary condition with topo- 
logical degree d E Z. Then, on the basis on conservation of vorticity, the limiting orbits for 
the gradient-flow dynamics show a vortex pattern. For q = 2, our main observation is the 
existence of vortices below the temperature of transition of phase for the mean fleld model, 
induced by the vorticity at the boundary of the lattice A, together with large ordered do- 
mains where the magnetizations ttt-a* become parallel. We discuss in detail dependence on the 
shape of the lattice, and on initial conditions. In particular, the application of the "simulated 
annealing process" allows the limiting configurations to move away from local minima, and 
reach lower energies. 

We also have some numerical evidence that, as in the case of Ginzburg-Landau functional, 
Kirchhoff-Onsager hamiltonian for the system of vortices gives a fairly good approximation 
of the minimizing free energy, despite of the non-local interactions. 

Finally, for q = 3, we examine in Sect. 5 the situation of spin- waves in the spirit of Belavin 
& Polyakov. 

Acknowledgements: We are very grateful to P. Piece who introduced us to the subject ; 
we also thank A. Messager and Y.Vignaud for many interesting and useful discussions. 
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1. Mean field approximation and renormalized Kac's Hamiltonian. 

Consider the lattice Z^, consisting in a bounded, connected domain A (the interior 
region), and its complement (the exterior region) A'^. In practice, we think of A as a large 
rectangle with sides parallel to the axis of Z^, of length of the form L = 2"^, n G N. Physical 
objects make sense in the thermodynamical limit A — > Z^, but in this paper we work in large, 
but finite domains. 

To each site z G Z^ is attached a classical spin variable ai G S*^"^, q = 2,3. The 
configuration space X{Z'^) = (S^~^)'^ is the set of all such classical states of spin ; it has the 
natural internal symmetry group 0~^{q) acting on S*^"^. The state a G A'(Z^) will denote the 
map cr : Z^ ^ S*^"^, i ^ cr(^)- Given the partition Z^ = A U A'^, we define by restriction the 
interior and exterior configuration spaces X{A) and X{A'^), and the restricted configurations 
by a\ and a\c. The Hamiltonian in Z^ describes the interaction between different sites 
through Kac's potential defined as follows. 

Let < J < I be a function on with compact support and normalized by J = 1- 

We can think of J also as a function on the lattice. There is a lot of freedom concerning the 

choice of J, but for numerical purposes, wc take J as f/2 the indicator function J of the unit 

rhombus with center at the origin, in other words J{x) = J(|x|i) where | • |i is the norm in 

R^. Thus the support of J is thought of as a chip of area 2, and considered as a function on 

the lattice, J takes the value 1 at the center, and 1/4 at each vertex, so that J{i) = 1. For 

ieA 

7 of the form 2""^, we set Jj{x) = 7^7(70;), and extend the definition above in the discrete 
case so that enjoys good scaling properties, namely the stratum of full dimension (i.e. the 
set of points interior to the chip) has weight 1, the strata of dimension 1 (the points on the 
sides on the chip) have weight 1/2, and those of dimension (the vertices of the chip) have 

weight 1/4. Thus, again ^ J-^ii) = 1- The discrete convolution on A is defined as usual. 

iez2 

For instance, {J^*a){i) = J-yii ~ represents, with conventions as above, the mean 

J6Z2 

value of cr over the chip of size 7"^ and center i, with a weight that depends on the stratum 
containing j. 

Note that we could replace the lattice Z^ by the torus (Z/LZ)^ or the cylinder (Z/LZ) x 
Z, which amounts to specify periodic boundary conditions in one or both directions. Ther- 
modynamic limit is obtained as L ^ cxd. 

The coupling between spin at site i and spin at site j is given by J-yii—j) ; this is known 
as Kac's potential. From Statistical Physics point of view, Kac's potential, for small 7, shares 
locally the main properties of the mean field, i.e. long range ~ 7~^, large connectivity ^ 7"^ 
of each site, small coupling constant a; 7^ of the bonds, and total strength of each site equal 
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to 1. 

Given the exterior configuration cta^ £ we define the Hamiltonian on as 

(1.1) H^{aA\(JAo) = ~^Mi-j){(j{i),<j{j))- M'^- j){^i^)^^U)) 

i,jeA {i,j)eAxA<= 

where cr{i), for simphcity, stands for a\{i) or crAc(*), and (•, •) is the standard scalar product 
in R^. We note that as J > 0, the interaction is ferromagnetic, i.e. energy decreases as spins 
align. 

We give here some heuristic derivation of the model we will consider, starting from 
principles of Statistical Physics. A thermodynamical system at equilibrium is described by 
Gibbs measure at inverse temperature We assume an a priori probability distribution v 
for the states of spin, and because of the internal continuous symmetry of A" (A), we take v 
as the normalized surface measure on S^""*^, i.e. v{dai) = u!~^S{\ai\ — l)dai, where Uq is the 
volume of S^"-*^. Then Gibbs measure on X{A) with prescribed boundary condition ctac is 
given by 

(1.2) /i/3,7(<i(TA|crAc) = —j^ -exp[-/3i/^(cTA|cTAc)] TT Zy((ic7A(i)) 

where Zp ^{a\c)^ the partition function, is a normalization factor which makes of |U/3,-y a 
probability measure on A'(A), conditioned by a\c e X{A^). It is obtained by integration of 
1^ 13, 'y {da aWa-) over Oq = (S*^"^)"^. 

Since we are working on Z^, there exists, for any /3 > 0, 7 > 0, an infinite volume Gibbs 
state fip^-y, i.e. a (unique) probability distribution fifs^-^ on the space X of all configurations ob- 
tained by taking the thermodynamic limit A — > Z^. This measure satisfies suitable coherence 
conditions, i.e. DLR equations. 

Nevertheless, we are faced with various difficulties, indicating that fJ- p^j [da aI^A'^) should 
not be the object to be directly considered. It is known that (and this goes back to Van Hove 
for the Ising ferromagnet, i.e. q = 1, see [Si, p. 31], ) in order to understand thermodynamical 
properties for spins models, one should instead average spins over mesoscopic regions and 
consider the image of Gibbs measure through this transformation, the so called "block-spin 
transformation" . So we introduce the empirical magnetization in the finite box A C Z^ 

(1.3) ^a(^) = ^E^W 

and given any m e R^, |m| < 1, we define the canonical partition function in A as 

(1.4) Zp^'^^m) = [ exp[-pH^{aA\(TAc)] TT iy{da(i))6{mA{a) - m) 

•^(s'-O ieA 
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see [Si,p.31]. Also for Kac's model, it can be shown, taking the thermodynamical limit 
A — > Z^, that the quantity 

(1.5) F^iP, m) = - limA^z^ ^ log Z^;^^" (m) 

is well defined, and doesn't depend on the boundary condition on ; it will be interpreted 
as the thermodynamic free energy density of the system. It is defined for a system with finite 
interaction of range 7"^, i.e. before taking the mean field limit 7—^0. 

So far, parameter 7 was kept small but constant ; the limit 7 ^ is called Lebowitz- 
Penrose limit. Let 

(1.6) /^(m) = -^|m|^ + -^/(m) 

be the free energy for the mean field, I{m) denotes the entropy, see (2.2) below. 
Lebowitz-Penrose theorem (in this simplified context) states that 

(1.7) lim^^o F^ip, m) = CE{fp{m)) 

See [BuPi] for the case of a 1-d lattice and continuous symmetry, the proof can be carried 
over to Z^. Here CE denotes the convex envelope, to account for Maxwell correction law. 

From this we sketch the renormalization procedure that leads to Lebowitz-Penrose the- 
orem, as stated e.g. in [Pr,Thm. 3.2.1] for q = 1, following earlier results by [AlBeCaPr] 
(actually, this is the "pressure" version of Lebowitz-Penrose theorem, but the argument can 
easily be adapted to free energy. ) The following result will not be used in the sequel, we just 
give it for completeness. 

Since we will take (in this paragraph) a continuous limit, we do assume that J is a 
differentiable function, not necessarily of compact support, but with ||VJ||i < 00 (the 
norm. ) The lattice dimension d can be arbitrary. We set, following (1.3), A = A(a:;) and 
^a{x) = ^X(x)^^-'' ^'^^'^ A{x) will be a square "centered" at a variable x G Z^, with sides 
of length ^, 6 of the form 2"^, p G N, ^ much smaller than the diameter of A, but still 
containing many sites, for instance diam (A(x)) = 7"^/^. Actually, we need to replace (1.3) 
by an integral, which allows to extend x 1— > maix) on R*^, but for simplicity, we present it as 
a discrete sum. 

The averages mcr{x) are called (empirical) magnetizations. The set of all such magneti- 
zations rrio- G R"^ is the image of A'(Z^) by the block-spin transformation yr-y : d — > m^, and 
will be denoted by A:'(Z^). This is the set of "coarsed-grained" configurations. 

It has again the continuous symmetry group 0^{q), and this is a subset of the convex 
set M. of all functions m : Z^ — > R*^ such that |m(a;)| < 1 for all x. When considering 
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microscopic interior and exterior regions as above, the partition = A U A*^ induces of 
course restricted configuration spaces X{A*) and X{A*^), where A* = {a; e : A(a;) C A} 
and A*'^ = {xeZ^ : A{x) C A^}. So let m e M. 

We introduce as in (1.4) the canonical Gibbs measure conditioned by the external con- 
figuration a Ac = cr*^ : 

(1.8) ii^^^,A{da;m\a^) = - — / Tl^{dcrii)) exp[-pH^{aii)\a^)]S{7rjaii) - m) 

where the partition function ^ a^'^a^^ defined in (1.2). For simplicity, we have removed 
the index A from a. By definition of the image of Gibbs measure through the block-spin 
transformation, we have 

(1.9) / dmiJ,p^j^Aida;m\a!y) = 1 

J\m\<l 

where dm is the normalized Lebesgue measure on the product space HaeA* ^qi^j i-^gi^j 1) 
denotes the unit ball of R'^. ) Let 



jF(m|m^) = - I dr I dr' J{r — r')\m{r) — m{r')\'^ 
4 -'Ao Jao 



(1.10) ^ 



+ - / dr I dr' J{r — r')\m{r) — m{r')f + j dr{^fi3{m{r)) — f(j{mi3)) 



' A^ 



be the continuous, free energy in a box Aq C of fixed size Lq, rescaled from A by a 
factor proportional to 7. Here mp is the critical value for the mean field see Sect. 2. 
Assume, as before, that the diameter of all block spins A(a;) equals 7"^/^. Then we can give 
a special meaning to the approximation /i/3,7,A ~ exp[— /37~''^(m|m'^)] (in the logarithmic 
sense) stated in the Introduction, by establishing the analogue of [AlBeCaPr, Lemma 3.2] in 
case of continuous symmetry, improving also [BuPi, Lemma 3.1]. Let e be any (fixed) unit 
vector in R^, and rhp the constant function on A equal to mpe, which we extend to be equal 
to rrf on A*^. We have the following : 

Proposition 1.1: Let q = 2. With the notations above, there are constants Ci, C2 > such 

that for any coarse-grained configuration m on A* : 

(1.11) 

-^(m)-(Lo7-')'(C2/3V7l|VJ||i + Ci7'/'log7-') 
< log[///3,7,A((icr; m|(7'')] -t- f3'y~'^J^{m\m'') 

< g{m) + T{mp\m'') + (Lo7-')''(C2/3V7ll V J||i + C^^'^'Hog^^) 

where g{m) = logna;GA*(l ~ \'m{x)\) 
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See [El-BoRo] for details. The divergence of g{m) as \m\ gets close to 1 reflects the 
fact that the entropy density I{m) is singular at |m| = 1, precisely where the mean field 
approximation breaks down, see also [BuPi, Theorem 2.2]. So the approximation /i/3,7,A ~ 
exp[— /57~^.F(m|m'^)] holds true when the magnetization stays bounded away from 1, as is 
the case in most applications. 

Having this construction in mind, we shall proceed the other way around, and make a 
simple renormalization of (see Proposition 2.1 below). Actually our sole purpose is to 
give a discrete analogue for the excess free energy functional as in (1.10), most adapted to 
numerical experiments on the lattice. 

2. Renormalized Hamiltonian on the lattice. 

We restrict here to g = 2, in Sect. 5 we show how these considerations easily extend 
to q = 3. Recall from (1.6) the free energy for the mean field, I{m) is the entropy function 
of the a priori measure u, which can be computed following [BuPi]. Namely, introduce the 
moment generating function 

(2.1) (j){h) - / e^'^'^'^duia) 
and define I{m) as Legendre transformation 

(2.2) I{m) = I{\m\) = supf,^j^,{{h,m) - log0(/i)) 

We denote by In the modified Bcsscl function of order n. For q — 2, we have (j){h) = (/>(|/i|) = 
Io{\h\). Function p I{p) is convex, strictly increasing on [0, 1], I{p) ~ as p — 0, 
/(p) ~ — ^log(l — p), as p — > 1, and these relations can be differentiated. We have also 
/' = ((log(^)') ^ and (log(/))'(t) = /i(t)//o(t), this is of course a real valued function. The 
phase transition of mean field type is given by the critical point of the free energy i.e. 
the positive root of equation /3mp = I^mp), which exists iff /3 > I"{0) — 2. So the critical 
manifold has again 0^(2) invariancc. 

Now wc specify the choice of mcsoscopic boxes A{x) and construct the excess free en- 
ergy functional by the procedure sketched above. When g = 2, it is convenient to use the 
underlying complex structure of A'(Z^), so we shall write (1.1), with obvious notations, as 

(2.3) /^^((JAlaA-) = ^ J7(i-i)cj(i)a(j)-Re ^ J-,{i - jW{iWU) 

(i,i)eAxA= 

We introduce in detail the mesoscopic ensemble averages, or coarse graining approximation 
to renormalize Hj. Let 5 > be small, but still much larger than 7, we take again S = 2~'^, 
for some m G N. We take for A(a;), a; e Z^, a square "centered" at x, of diameter ^, 
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and of the form As{x) = {i = (11,12) E Ti'^ : ik ^ Tfi^ki^k + 1[} where we define as in 
(1.2), 777,5(0;) = (^)^ ^ cr(z). Thus we magnify by a factor 5/7 the "coarse graining" (or 

mesoscopic ensemble) labelled by a; e Z^, to the "smooth graining" (or microscopic ensemble) 
labelled by z e Z^. We have : 

Proposition 2.1: There is < a < ^ such that 

(2.4) (^fH^{aA\aAo) + U^^,{ms) + t/int(m5) - Wp{mp) = T{ms\ml) + |A|C»(52«) 



where 



(2.5) 



T{m5\mf) = ^ J5{x -y)\ms{x) -ms{y)\^ 

x,yeA* 

+ ^ Js{x-y)\ms{x) -ms{y)\'^ + Y fp{ms{x)) - ffiimp) 

{x,y)eA* xA*'' xeA* 

Ue^tims)^^ ^ Jsix -y)\ms{y)\'^ 

(a;,j/)eA*xA*'= 

Uintims) = ^ Y I{m5{x)) 



Proof. To start with, consider the first term in (1.1) 

"Y - j)(^i^)(^U) = Yl J5{x-y)m5{x)m5{y) 

i,j&A x,yeA* 

+ 7' E Y^ {J{^{i-3))-J{d{x^y)))a{^)W) 

and denote by i2(A*) the second sum in the RHS of (2.6). Let Cq = Bi(0, J) be the rhombus 
(or £^-ball in R^) of center and radius |, corresponding to the shape of the interaction J, 
and for e Z^, its translate C^' = ^x' + Co, we denote also by C*, C A* the corresponding 
lattice obtained from Cx' by deleting 2 of its sides, so that A* = U^'ez^ ^x' (disjoint union), 
and A* is covered by those C*, with x' = {x'l, x'2), x'j G {±1, • • • , ±7-£'}- Let also E{x, y) = 
e As{x) X As{y) : J{'y{i — j)) — J{S{x — y)) 7^ 0}. We can consider E{x,y) as a 
symmetric relation E : A* ^ A*, E{x) = {y G A* : E{x,y) ^ 0}. By translation invariance 
of J, for any x' e Z^, we have \E{x,y) \ = \E{x — \x' ,y — jx')\, so that 

(2.7) \E{x,y)\<A(ff Y 

x,yeA* x,yeC* 
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With the choice of £^ norm, we have E{x,y) ^ for all x,y e Cg, and maxxec* = 
(l + ^)^, while mmx£c* \E{x)\ is of order unity. In any case, |-E'(a;)| depends on x and 5, 
but not on 7, and it is easy to see that for some < a < = 0{6~'^^^~°'^), 

5^0. [Actually, this kind of estimate is well-known, see e.g. [BlLe] and references therein 
for related results, and applies whenever the support of J is a convex set. ] 

On the other hand, we have the rough estimate \E{x^y) \ < \A§{x) x As{y) \ = (^)^, and 
since \(7{i)\ = 1, 

I E J2 ^ {J{^{i-3))-J{d{x-y)))a{i)W)\ 

^ xec* ^ 

This, together with (2.7), shows that i?(A*) < Const. (5^"L^. A similar argument gives an 
estimate on the remainder i?(A*|A*'^) for the second term in (1.1). Once we have replaced 
(5)^ E by E "^^^^ ~ y)'ms{x)ms{y) modulo R(A*) and i?(A*|A*''), which 

verify the estimate given in (2.4), we use the identity 

-2Rems{x)ms{y) = \ms{x) - ms{y)\'^ - \ms{x)\'^ - Imsiy)]"^ 

and express the "density" term ||toP in term of the mean field free energy fj3{m) as in (2.1). 
Summing over {x^y) and making use of the fact that J5 is normalized in l^{7?) eventually 
gives the Proposition. X 

Remarks: 1) In homogenization problems, one usually associates the discrete configuration 
a G X{K) with the function on taking the constant value a{i) on the square "centered" 
at 7^, i = {iiji2), i-e. on [7^1,7(^1 + 1) [x [7^2, 7(^2 + !)[• Furthermore the size of the domain 
A is normalized, so that taking the thermodynamic limit A — 00 is a problem of convergence 
for piecewise constant functions (or discrete measures) in some suitable functional space. 
As wc have seen in Sect.l, it is convenient to take a smooth interaction J. Thus a version 
of Proposition 2.1 was obtained in [BuPi] by replacing the discrete average ms{x) around 
X G A by an integral, or in [DeMOrPrTr] , [DeM] , [Pr] , ... by averaging over boxes of type 
Cx' as above. (For short we refer henceforth to the review article [Pr]). Since our ultimate 
purpose here consists in numerical simulations on a lattice, we chose instead to give a discrete 
renormalization for H^. 

2) Our renormalized Hamiltonian is now given by J-'{'ms\mg), the quantities we have 
subtracted are —Uex.tinT's) , attached to the configuration space A'(A'^), and —Uintims) that 
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can be interpreted as times the entropy of the system in A. Note we have also included 
self-energy terms i = j in the original Hamiltonian. Of course, relevance of this free energy 
to Gibbs measure (or rather its image through the block-spin transformation) after taking 
the thermodynamic limit, is a rather subtle question which will not be discussed here, since 
we content to finite lattices. 

3. Euler-Lagrange equations and non local dynamics. 

We are interested in the critical points of J-'{ms\mg). Denote as usual resp. by dm and 
dm the holomorphic and anti-holomorphic derivatives, we have for m = ms (for short), and 
any tangent vector of type (1,0) in the holomorphic sense. 



{dm^iTn\m'^) , dm) — ^ Js{x — y)(m{x) — m{y))Sm{x) 

{x,y)eA* xA*^ 

+ ^ XI Mx - y){rnix) -m(y))5m{x) + ^ {~m(x) + ^ ^^^^\ x))5m{x) 



x,y&A* xGA* 

Using again the normalization of J5 in £^(Z^), the relation I{m) = /(|m|), and setting as 
before Js * m{x) = Js{x — y)m{y), we obtain 

(3.1) {dmJ^{m\m''),Sm) = I Yl {-Js*m{x) + ^ [^^^"^ m{x))dm{x) 
Since is real, this gives Euler-Lagrange equation : 

(3.2) -Js * m{x) + \ ^-^^^^m{x) = 

p \m\ 

Let / = {J-')~^ = ^ denote the inverse of the function Thus / : [0, -|-oo[^ [0)1[ is 
strictly concave, /(O) = 0, /'(O) = 1/2, and /(p) ^ 1 as p — >^ -\-oo. Since the inverse of 
m /'(|m|)|^ defined on the unit disk is given by n 1— > /(l^l)]^, n e C, (3.2) takes the 
form 

(3.3) -m + /(/3| * m|)f^^ = 

Following [Pr], to find the critical points minimizing the excess free energy functional T we 
solve the "heat equation" 

,3.4) =|i = _ + ;(^|,,,„,|)^i„A. 
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with prescribed (time independent) boundary condition on A**^, and initial condition m\^* = 
mo- By Cauchy-Lipschitz theorem, equation (3.4) has a unique solution, defined for all t > 0, 
valued in X{A*). Monotonicity of is given in the following : 

Proposition 3.1: is a Lyapunov function for equation (3.4), i.e. there exists a free energy 
dissipation rate function X : A" (A*) — > R"*", X(m) = iff m solves (3.3), and 

±J^{m{;t)\m-)=-l{m{;t)) 

along the integral curves of (3.4). 

Proof. We have, using (3.1) and (3.4) 
(3.5) 

X(m(-,t)) = -— = -{dmJ^, —) - (BmT, —) 

= \^^Y. i-f^-^s * ^(^) + ^-^TT^Mx)) {m{x) - f{(3\Js * m\) f/' * . (x)) 

Let m = pe*^, /JJj * m = p'e*^', T(rn{-, t)) equals a sum of terms of the form 

R = |(p7(p') + pT'ip) - {PP' + f{p')i\p)) cos{o - e')) 

then using (p — f{p')) {l'{p) ~ p') > for any p, p' since /' is increasing, we obtain the lower 
bound > 1(1 - cos(^ - e')) {pp' + f{p')T'{p)) > 0. And because pp' + f{p')I'{p) = iff 
p = or p' = 0, this estimate easily implies the Proposition. X 

From Proposition 3.1 and a compactness argument as in [Pr], follow that in the closure 
of each orbit of equation (3.4) there is a solution of (3.3), or equivalently, of Euler-Lagrange 
equation (3.2), i.e. a critical point for T. As suggested by numerical simulations, this 
critical point is not unique, and depends on initial conditions (except of course when j3 < 2.) 
We expect however some uniqueness in the thermodynamical limit A* — > oo, modulo the 
symmetry group. 

Now we give estimates on solutions of (3.4) or (3.3), borrowing some ideas to [Pr]. Eq. 
(3.4) can be rewritten in the integrated form : 

Tx * 771 

(3.6) m(x,t)^e-^m(x,0)+ dtie^''^ fipUg * m\)-^ i{x,ti) 

Jo iJs*^] 

An effective construction of the solution is given by the "time-delayed" approximations. It 
will also be used, discretizing time, in the numerical simulations below. We define inductively 
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mh{x,t), /j. > 0, on the intervals [hk,h{k+ 1)[, A; e N, by mh{x,t) = e ^mo{x) ior < t < h, 
and for hk < t < h{k + 1), k > 1 : 

(3.7) mh{x,t) = e^^-*mh{x,kh) + f dhe*'-* fiPe'^lJg * mh\) * , {x,ti - h) 

Jhk \Js*mh\ 

Using Lipschitz properties of the coefficients, it is easy to see that, as /i — 0, mh{x,t) tends 
to the solution m.{x,t) of (3.4) uniformly for x E A* and t in compact sets of R-|-. We prove 
estimates on m{x,t) using sub- and supersolutions of (3.4). We start with : 

Lemma 3.2: Assume /? > 2, and let A(t), t > be the solution of 

(3.8) ^{t) + X{t) - f{pX{t)) = 0, A(0) = A e [0, 1[ 
If A > m/3, then X{t) < A for aU t > 0. 

Proof. Write (3.8) in the integrated form as in (3.6) and consider the approximating sequence 
Xh{t). Since Xh{t) tends to X{t) uniformly on compact sets of R-|-, it suffices to show the 
property stated in the Lemma for A/i, and h > small enough. For < t < h, Xh(t) — Xe~*, 
so the property holds, while for h < t < 2h, performing the integration in (3.7), we get 
Xh{t) = e-*A + Jldtie^'-^f{/3e-*'X). Since x > f{Px) iff Px < T'{x) (whence iff x > mp, ) 
if A > m^, and h > small enough, then /(/3e~*iA) < e~^^X, and Xh{t) < X. By induction, 
using also that / is increasing, but without changing h > anymore, it is easy to see that this 
property carries over for all t > 0. By a continuity argument, this holds true for all A > m^. 
«■ 

Then we claim that the modulus of the magnetization doesn't increase beyond mp. More 
precisely we have : 

Proposition 3.3: Assume /3 > 2, and let m{x,t) be the solution of (3.4) such that mo{x) = 
m{x, 0) satisfies |mo(a;)| < A < 1, for some A > mp, and all a; G (so including the boundary 
condition on the exterior region .) Then \m{x,t)\ < X for all a; e Z^, and alH > 0. 

Proof: Eq. (3.6) shows that 

(3.9) \m{x,t)\ <e~^\mo{x)\+ [ dtie^'~^ f{(3\Js * m\){x,ti) 

Jo 

Now by the monotony properties of the convolution and the function /, we have f{P\Js * 
m|)(x,ti) < f{PJs * |w|)(a;,ti), so the solution X{t) of (3.8) with A(0) = A is a supersolution 
for (3.9), and Lemma 3.2 easily implies the Proposition. X 

We now look for lower bounds on m{x, t). Since there are in general vortices, one cannot 
expect a global, positive lower bound on \m{x,t)\, unless there is no vorticity on initial and 
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boundary values. On the other hand, we know (at least for a 1-d lattice, see [BuPi] , ) that the 
Gibbs measure of the configurations at equilibrium ms G A'(A*|A*'^) with |m5(a;)| arbitrarily 
close to mj3, has to be large. We have : 

Proposition 3.4: Assume > 2, and let m{x,t) be the solution of (3.4) such that mQ{x) — 
m{x, 0) as in Proposition 3.3 satisfies Re(z/mo(a:)) > /U, for some fixed z/ G ~ {2; G C, \z\ — 
1} and fi > and all x G Z^. Assume furthermore that fi satisfies (fi"^ + A^)^/^ < l3f{j3X), 
where A is as in Proposition 3.3. Then Re(z/m(x,t)) > ^ for all x G A*, and all t > 0. 

Proof. As in the proof of Proposition 3.3 we make use of a comparison function. So let jiit) 
verify the differential equation 

(3-10) f (0+MO-/3/(/3A) ^^^^^,^|^^^^^ M0)>/. 

Write (3.10) in the integrated form as in (3.6) and consider the approximating sequence 
jJih{t) as in (3.7). We shall show that jihit) > A* for alH > provided /i(0) > ji verifies the 
inequality given in the Proposition. Namely, this holds for < t < /i, because then iih{t) = 
//(0)e~* > fi ioT h > small enough, while for h < t < 2h, performing the integration as in 
(3.7), we get iJ,hit) = e'^-V(O) +/?/(/?A) dtie*^-^ iie'^^ {{iie-^^f + A2)-V2. By hypothesis, 
fie~*^ Pf{PX){(fie~^^)^ + A^)~-^/^ > for h small enough. So again fih{t) > M- By induction, 
using that the function p 1— > p(p^+A^)~^/^ is increasing on R+, it is easy to see that iJ,h{t) > jJi 
holds for all t > 0. Because the coefficients of (3.10) are uniformly Lipschitz, jJih{t) tends to 
p,{t) uniformly on compact sets in R+, and this property holds again for /i(t). 

Now we turn to the equation for m{x,t). Possibly after rotating the coordinates, we 
may assume 1^=1, i.e. Re mo (a;) > fx and all a; G A* (again, we have included the boundary 
condition in the initial configuration. ) Write m{x,t) = u{x,t) + iv{x,t), u,v real and take 
real part of (3.4). The integrating form of the resulting equation writes : 

(3.11) uix, t) = e-'uix, 0) + r dhe'^-'f{(3\Js * mh ^''. (x, h) 

Jo (3\Js*m\ 

As p' I— > is decreasing on R+, and by Proposition 3.2, [Jg * m| < (jJs * itp + A^)"*^^^, we 
have 

mJs * m|) ^ f{p{\Js*u\^ + X^f') ^ fipX) 



p\Js*m\ /?(|J5*«|2 + A2)'/' /3(|J5*«|2 + A2)'/' 

the last inequality because / is increasing. Since u{x, 0) > p,, by continuity we have u{x, t) > 
at least for small t > 0, and (3.11) gives 

(3.12) u{x,t)>e-^u{x,0) + Pf{(3X) I dtie*'-\Js * u){{Js * uf + X'^y^^^{x,ti) 

Jo 
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Now, using the monotony of the convolution, and again the fact that the function p i— > 
p{p^ + A^)"-*^/^ is increasing on R+, we can easily show that the solution ji{t) of (3.10) with 
//(O) = jiis actually a subsolution for (3.12), for alH > ; the properties proved already for 
IJL{t) then imply the Proposition, ijt 

Of course, by continuity. Propositions 3.3 and 3.4 imply the corresponding estimates for 
the solutions of (3.3), or equivalently for the solutions of Euler-Lagrange equation (3.2). Our 
last result states that if /3 < 2, then m(x, t) tends to t — oo, which is consistent with the 
absence of phase transition (or spontaneous magnetization) at high temperature. 

Proposition 3.5: Assume /? < 2, and let m{x,t) be the solution of (3.4). Then m{x,t) 
on A* as t ^ +oo. 

Proof. Using that f{p') < \p' , all p' > 0, (3.9) shows that 

\m{x,t)\ <e~^\mo{x)\ + - / dtie^'-'^Js * \m\{x,ti) 

2 Jo 

So by taking convolution 

3 f^^ 

Js*\m\{x,h) <e-^'Js*\mo\{x) + ^ / dt2e*^"^Vf * |m|(a;,t2) 

2 Jo 

and integrating the resulting inequality : 

\m{x,t)\ < e-* [|mo(x)| + ^Js * \mo\{x) + (^) V^) (e^') Jf * \m\{x, ■)){t)] 

where T^^^u{t) = dti J^^ dt2 ■ ■ ■ dtku{tk) denotes the /c-fold integral of u, and J^'^ the 
/c-fold convolution product of with itself. By induction, we get : 

g 3 kt^ 

\m{x,t)\ < e-*[|mo(a;)| + ^tJs * \mo\{x) + ■ ■ ■ + (|) -j^J^ * \mo\{x) 

+ r('=+i)(e(-)j;('=+'^*|m|(x,-))(t)] 
The series is uniformly convergent for t in compact sets so we can write 

+°o P ,fk 

|m(x,t)|<e-*J](|)'-jf *KI(a:) 

k=0 

When /3 < 2, using Jg * \m\{x,0) < |?Bo(ic)| ^ 1, it follows that m{x,t) — > for all x G A* 
as t — s> oo. This holds again for f3 = 2 since we may assume that mo has compact support, 
and we know (see [Ho, Lemma 1.3.6]) that J^^ uniformly on (or on in the discrete 
case, ) as A; — >• oo. X 
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4. Vortices. 

We consider here the problem of finding numerically the critical points of Euler-Lagrange 
equation (3.3) by solving (3.4) subject to a boundary condition on A**^ presenting vorticity. 

a) Generalities. 

First we recall some facts about the degree of a map. Let m : — > C be a differentiable 
function, considered as a vector field on R^, and subject to the condition |m(a;)| — > £ > as 
|a;| — > oo uniformly m.x = x/\x\. Then the integer 

(4.1) deg^m=;^/" d(argm) = / — 

J\x\=R ^™ J\x\=R ^ 

is independent of R when i? > is large enough, is called the (topological) degree of m at 
infinity, and denoted by deg^^ m. 

We define in the same way the local degree (or topological defect) deg^^ m oi m near 
xq, provided m{x) x ^ xq, by integrating on a small loop around xq. The local degree 
takes values dj e Z. When m has finitely many zeros Xj inside the disc of radius i?, its 
total degree (or vorticity) is defined again as the sum of all local degrees near the x^'s. In 
many boundary value problems (or generalized boundary value problems, in the sense that 
the boundary is at infinity, ) such as Ginzburg-Landau equations, total vorticity is conserved, 
i.e. degg^m = ^ deg^^ m. Generically dj = ±1 ("simple poles".) Our aim is to check this 
j 

conservation principle in the present situation. 

We can define analogously the degree of a discrete map, which makes sense at least in 
the thermodynamical limit. If m{x) — p(x)e'^^^^\ the degree of m at infinity is the degree 
restricted to the lattice A*^, e.g. by 

(4.2) d = degA*c J2{0j+i - Oj) 

j 

along some closed loop Fj, C A**^ encircling A*, the sites along Fj, being labelled by j, assuming 
that this integer takes the same value on each Fj,. 

The local degree near xq, where m{xQ) = 0, is identified again by computing the angle 
circulation on a loop encircling xq. Local degrees are also expected to take, generically, values 
±1. 

We chose our parameters as follows. We start with prescribing the degree of the spin 
variable a on A*^, and take on F,,, the i:th loop away from A, containing sites, (A^^, = At+P, 
where P is the perimeter of A, we take enough t's to cover the range of interaction, ) with a 
uniform distribution: 

(4.3) aj = exp t{2TTdj/N, + 0o), 1 < J < iV. 



16 



here 0o is a constant (e.g. 0o = 1) that "breaks" the symmetry of the rectangle A. We shall 
also randomize these boundary conditions. 

To this spin distribution on A'^, we apply the block spin transformation (1.3), so to have 
a distribution of magnetization on A*^, then we prescribe initial conditions inside A*. The 
simplest way is to take zero initial values, which gives a particular symmetry to the solution. 
Otherwise, we can choose them as random numbers, either small, or with absolute value less 
than rri/j. All these cases will be discussed. 

We usually fix the inverse temperature /? = 5, so = 0.72 ; the results do not depend 
on P in an essential way, we just observe that magnetization tends to as /? — > 2+. The 
diameter L of the lattice A ranges from 2^ to 2^°, the size 5/^ of the diameter of the block- 
spin A (a;) is set to 4 (most of the time) so the diameter L* of the lattice A* ranges from 2^ 
to 2^. The lattice is either a square, or a rectangle. 

The size 1/5 of the length of interaction in A* ranges from 2 to 32, thus the corresponding 
interaction in A has length I/7 = 4/d between 8 to 128. 

Equation (3.4) is solved by "time-delayed" approximations as in (3.7), implemented by 
the second order trapezoidal method to compute the integrals. 

These experiments lead to the following observations, vortices display in a different 
way, according to the initial configuration on A*, but always obey the conservation of total 
vorticity. 

b) Some typical configurations. 

We consider here the case of a uniform distribution of spins on the boundary. 

The particular case of zero initial values and a square lattice, gives raise to interesting 
symmetries (or degeneracies) in the picture : namely, vortices tend to occupy most of A* so 
to cope with the symmetry of the square. So for d = 1 there is a single vortex in the center, 
for d = 2 {cf Fig l.a) a vortex of multiplicity 2, (unless the degeneracy is lifted and turns 
into 2 nearby vortices, ) for d = 4 — 1, (cf Fig 2. a) one vortex of degree -1 surrounded by 4 
vortices of degree +1 near the corners, for d = 4, 4 vortices of degree +1 near the corners, 
for d = 4 + 1, same configuration as for d = 3, for d = 4 + 2 the picture looks alike, with 
a double vortex at the center, for d = 2x4 — 1, 4 new vortices appear near the center (cf 
Fig 3. a), etc. . .So the configuration depends essentially of the residue of d modulo 4 : new 
vortices show up from the middle towards the corners along the diagonals of A*. 
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Fig l.a: L* — 128, d — 2, zero initial condition 



Fig l.b: L* — 128, d — 2, random initial condition 




5 ID 15 3] a 3D 35 '-5 5 ID 15 3] 25 3D 



Fig 2. a: L* — 128, d — 3, zero initial condition Fig 2.b: L* — 128, d — 3, random initial condition 

Next we consider the case of a square lattice, but with random initial conditions, that 
is, we pick initial magnetizations with random direction and random length, provided the 
length is much smaller than mp, typically |mo(a:)| < 0.05. Then vortices are simple (i.e. have 
local degree ±1, total vorticity is of course conserved, ) and tend to display at the periphery 
of A* , in a pretty regular way, leaving some large ordered domain near the center. 

Thus, these configurations maximize the area of the lattice where the magnetizations are 
aligned , with an absolute value close to m^, (in accordance with the fact that energy 
decreases as the spins align. ) Their direction, in general, points out along one of the diagonals 
of A*. This is illustrated in Fig. 1,2, 3. b above, for a vorticity d = 2,3,7 respectively. In 
particular. Fig. 2 shows the topological bifurcation from d=4-l to d=3. These simulations also 
suggest that the equilibrium configurations depend on the initial conditions, but exceptional 
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configurations due to symmetry, for zero initial conditions, are essentially removed as soon 
small disorder is introduced. 
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Fig 3. a: L* — 128, d — 7, zero initial conditio 



Fig 3.b: L* — 128, d — 7, random initial condition 



Now we vary the shape of the lattice, changing the square into a rectangle, keeping in 
mind that thermodynamic limit, most of the time, should be taken in the sense of Fisher, 
i.e. the length of the rectangle A* doesn't exceed a constant times |A*|^/^. As expected, 
vortices tend to align along the largest dimension, but again, limiting configurations depend 
on whether the initial condition inside A* is set to zero or not. 
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Fig 4.c: Zy* — 256, i* — 128, d — 3, zero initial condition 

Thus, for zero initial condition, vortices display along the largest median of A*, with 
possible extra vortices near the corners (inheriting the features of the square lattice. ) Namely, 
they tend to repel each other so the energy cost in clustering is minimized by occupying the 
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corners. Typically, such configurations occur if d > 4 and the length of A* is only twice its 
width. But for sufficiently long lattices, or small degree, they just stand the median line. See 
Fig. 4. a and 4.c. 

For small random initial conditions as above (Fig.4.b), we recover the general picture of 
square lattices, i.e. vortices set along the boarder of A*, leaving a large space in the middle 
with parallel magnetizations. In any case, degeneracies are lifted, and all vortices have degree 
+1. 

c) The simulated annealing. 

If we increase the initial conditions, still keeping |mo(a;)| < mp, we obtain similar pic- 
tures, but with a non uniform distribution of defects : conservation of total degree holds, but 
at the same time, many vortices spread over the lattice, and the corresponding long range 
order region shrinks correspondingly. This suggest that the gradient-flow dynamics converges 
only to a local minimum of the free energy. 

For reaching lower energies, we let the system explore other regions of the configuration 
space. This can be achieved through simulated annealing, see e.g. [KiGeVec]. Replace the 
dynamics (3.4) by 

(4.4) ^ = + fm\Js * in A* 

at \Js*m\ 

where /3(t) depends continuously on t G [0,ti], starting with (3o < f3i = /3, with negative 
slope at t = 0, so that the system is heated initially up to a peak ^ 1/2 (the critical 
temperature) around t = t2, and then gradually cooled down to (3 aX t = ti. Function (3{t) is 
oscillating between successive warm and cool periods, so to "shake" sufficiently the system. 
Then we keep the temperature constant till we reach equilibrium. 

It is not difficult to optimize, empirically, the annealing function /9(t), and our choice 
was the following : 




Fig 5: the annealing function 
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We applied this method first to the case of a square lattice, when the equilibrium config- 
uration corresponding to some total degree do is used as an initial condition for a dynamics 
with degree di. We fix Pi = 5, L* = 128. 



Consider first the case di = 3, the equilibrium configuration, with initial condition, is 
given in Fig. 2a, and the corresponding free energy is E = 99. We use simulated annealing 
to compute the equilibrium, starting from do = 4, —3, 5, and find respectively E = 23, 53, 51, 
see Fig. 6. So the energies obtained this way are less than with zero initial conditions, though 
the initial magnetizations are rather large. At the same time, symmetries get lost. Thus the 
cost for the 3 vortices to be aligned along one of the diagonals of the square as in Fig 6.b is 
less than to form a domino near the center as in Fig. 2a. 



35 

m 
W 

m 

Wi 
'Mi 

-5 




'''''''' ^ 'iir^I^-^-^N^^ 




10 



15 



20 



25 



30 



35 




10 



15 



20 



25 



30 



35 



Fig 6.a: rfo = 4, rfi = 3, = 23 



Fig 6.b: do = —3, (ii = 3, i? = 53 
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Fig 6.c: do = 5, (il = 3, i? = 51 



Fig 7. a: d\ — 5, (i — 7, zero initial condition 
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Fig 7.b: (Iq — 3, (ii — 5, without annealing 



Fig 7.c: (ig — 3, (ii — 5, with annealing 



In Fig. 7 we show how to pass from do = ^ to di = 5. The configuration with zero 
initial condition and di = 5 is given in Fig. 5a, and energy is i? = 113. Taking instead 
the equihbrium configuration for do = 3 as an initial condition yields, without simulated 
annealing, to Fig. 7b, with di = 6 — 1, and E = 216. Using simulated annealing gives instead 
Fig. 7c, which looks like Fig. 7a, and corresponding energy E = 115. Actually, the 3 vortices 
on the anti-diagonal of the square in Fig. 7b collapse into a single one at the center. 

Note also that the degeneracy in case of (ii = 2 (a vortex of multiplicity 2 at the center 
for zero initial condition, E = 24, ) is lifted through annealing from do = 4 : the 2 vortices 
move far apart, and E = —19. Other applications of simulated annealing will be given in the 
next subsection. 

d) More general configurations. 

We examine here the role of random fiuctuations in the distribution of spins on the 
boundary A*'^, so to account for possible defects in the structure. With notations of Sect. 4a, 
we take aj = exp(2iTTd{j /N^ + where e^j are uniform i.i.d. random variables with 

= 0, and (£(,,j)(,,i<j<Ar,-i, and variance small enough. The total degree is still equal 

to (i, but the variation of the direction of spins at the boundary is not uniform. As expected, 
the picture does not depart drastically from the previous cases. Vortices change their place 
according to the initial value, and tend again to gather inside A*, but take always the value 
+1 (assuming d > 0.) The sole effect of randomness in the boundary condition is to change 
the place of the vortices: namely they tend to get even closer to the boundary, so to leave 
larger ordered regions in the middle. 

In Fig 8.a,b, we have shown equilibrium configurations, obtained for d = 7, from the 
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same initial and boundary conditions, but with (resp. without) simulated annealing. Initial 
magnetization has been chosen at random, but a priori larger than before, the sole requirement 
being that |mo(x)| < m^. Random fluctuations on the boundary have been prescribed as 
above. 
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Fig 8. a: d, — 8 — E — 277, without annealing Fig 8.b: d — 7, — 240, with anneahng 

e) The KirchofF-Onsager correction. 

Another interesting result concerns the value of energy for the minimizing conflgurations. 
In case of Ginzburg-Landau equation, —Aif; + — = 0, where i{j is subject to a 
boundary condition with vorticity, it is known that energy of the minimizer vs. vorticity, has 
an asymptotic, as the n vortices Xj become distant from each other, the leading order term is 

n 

given by a "proper energy" , proportional to d'^ , and the next correction is the inter- vortex 

i=l 

energy given by so-called Kirchhoff-Onsager hamiltonian, of the form 
(4.5) Wq = —TV didj log \xi — Xj \ 

(see e.g. [BetBrHe] and [OvSi2] for precise statements. ) It can be interpreted as the electro- 
static energy for a system of charges dj interacting through Coulomb forces. It turns out that, 
despite forces in action have no electrostatic character, Kirchhoff-Onsager correction holds 
with a good accuracy in our case, even for long range interactions (i.e. for small 7, ) but 
provided the inter-vortex distance is bounded below by the range of the interaction. We have 
listed below some graphs of X = jF(-|m'^) — Wq, obtained with uniform boundary conditions, 
which show that K roughly grows linearly with d (cf Fig 9). 
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Fig 9. a: L* — 128, zero initial condition 



Fig 9.b: L* — 128, 3 random initial conditions 



Fig 9.b shows that several random trials for initial conditions give approximately the 
same renormalized energy K. 

5. The Heisenberg model 

We consider here "stationary spin waves" for = 3, in a setting similar to this of Belavin 
& Polyakov [BePo], [Po,Chap.6]. 

Let us first recall the model. We look for minimizers of H{a) = \'Va{x)\'^dx, among 
all configurations a : subject to the condition a{x) (0,0, 1) as \x\ oo. This 

boundary condition not only ensures a finite energy on the whole plane, but also allows to 
extend a as a map on the one point compactification of R^, so we may consider its degree 
D{a) G Z, or winding number, on the sphere. Differentiable maps are classified by 

their degree, in the sense that cr, a : — are homotopic iff they have the same degree. The 
main result of Belavin and Polyakov asserts that there exist solutions of that minimization 
problem, with given degree ; they are called instantons, and expressed (in suitable coordinates 
associated with the stereographic projection C — > S^) by arbitrary meromorphic functions of 



the form Yl 



d 



Here (oj, hj) E C play the role of vortices in the case q 



2 ; they have 



i.j = l z-bj 

a natural structure of dipoles, with poles placed at aj and bj. So the minimization problem 
(for a given homotopy class) has a continuous degeneracy, parametrized by the family (a^, bj) 
which we interprete as moduli. The energy of all such instantons is a constant proportional 
to D. 

It is then natural to consider the contribution of all instantons of same energy D. Some- 
what heuristically, [Po] obtains, after summing over D G N, a grand partition function of the 
form 

(5.1) E:(A) = ^ -— ^ / Yldajdbjexp[^{log\ai-aj\^ + log\bi-bj\^)-^log\ai-bj\'^] 

D>0 ^ 



and each instanton behaves as if it consisted of a pair of opposite Coulomb charges, placed 
at ttj and bj. Since the 2 dimensional Coulomb energy is given by (l/47r) log|aj — the 
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exponent in (5.1) reminds us of the Kirchoff-Onsager hamiltonian (4.5), and formally, S(A) 
is the grand partition function of a plasma at inverse temperature (3 = Att. 

It is not known to which extend these instantons are stable relatively to perturbations 
of H{a), e.g. due to the influence of temperature. 

We start with some considerations on the degree of a map on Z^. Let m : — > 
be a discrete map, defined through the stereographic projection — > S^, the one point 
compactification of Z^ given by Z f« Z^ U {oj}. The coordinates on the source and target 
space are given by the polar and azimuthal angles {d,(fi), and {6,(f) respectively. 

Consider the complex C = (Z^, Lz2, Pza) and its homology group. Here Lz2 is the set 
of bonds of unit length indexed by closest neighbors x,x' e Z^, and Pz^ the set of chips of 
unit area (plaquettes) around a; e Z^. See e.g. [A] for concepts of polyedral topology. 

We define as usual the discrete jacobian Jac m(a;) = computed on the plaquette 

around x. Let yo = m{xo),xo e Z^ be a regular value of m, i.e. Jacm(a;o) ^ 0. The integer 

D^g{m) = \{x e m~^{{yo}) : detJac m(a;) > 0}| — |{a; e m~^{{yo}) : detJac m(a;) < 0}| 

is called local degree of m at xq. In case where Dj^^^im) takes the same value for all Xq G Z^, 
we call it the degree of m and denote by D{m). This is the general case, and D{m) counts 
the number of coverings of the sphere. Then D{m) will be given by the discrete analogue of 
the integral 



computed on the complex C defined above. When the values of m avoid a neighborhood of 
cd, we put D{m) = 0. If Dx^f{m) = d for all Xq in a neighborhood of G Z^, we call d the 
degree of m at infinity and denote d = Di^{m). See e.g. [BlGaRuSh] and references therein 
for a more complete study of topological properties of discrete maps. 

We conjecture that for Kac-Heisenberg model, if m is a minimizer for the free energy 
J-'{-, \m'^), i.e. m solves (3.2) or (3.3) with D^{mo) = d, after we take the thermodynamical 
limit A — s> Z^, then either m vanishes at some point a; G Z^, or : Z^ — > has degree 
D. In practice however, we have only observed configurations with < D < d. So m 
shares some features with Belavin & Polyakov instantons, though with less symmetries or 
degeneracies, and a possible "degree loss" from infinity, since we are not really working in 
the thermodynamical limit. 

It is straightforward to extend the constructions of Sect. 2 and 3 to the case q = 3. Let 
us sketch the main steps. The moment generating function is now 4'{h) = (p{\h\) = '^^"^| , see 

[BuPi], and for the entropy function I{m) = I{\m.\) defined in (2.2), we have/' = {{log (/))') , 
and (log0)'(t) = L(t) = — j (the function / before) is known as Langevin function. 
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This is a concave, increasing function on R"*", L(t) ~ t/3 as t ^ 0, and L{t) — 1 as t — oo. 
There is a phase transition of mean field type i.e. a positive root for equation (3mp = I'{mp), 
iff /9 > /"(O) = 3. We derive Euler-Lagrange equations for J^{ms\mf), as in Sect. 3 (here we 
simply see m as a vector in the unit ball of R^, the complex representation of m was not 
essential, ) and find 



(5.2) 



-m + L{(3\Js * m\)- 







\Js * m\ 

For the corresponding gradient-flow dynamics (3.4), there is again a free energy dissipation 
rate function, which we compute exactly as in Proposition 3.1. Furthermore, we have esti- 
mates on m{x,t) as in Propositions 3.3 and 3.4 ; more precisely 

Proposition 5.1: Assume /9 > 3, and let m{x,t) be the solution of (5.2) such that mo(a:) = 
m(a:,0) satisfies |mo(a;)| < A < 1, for some A > m^, and all x G Z^. Then |m(x,t)| < A 
for all X G A*, and all t > 0. Assume moreover the z-component ttiq^x) of mo(x) satisfies 
Wo(x)) > ^ > 0, for aU x G Z^, and some // > with (//^ + A^)^/^ < PL{j3X). Then 
m''(x, t)) > /i for all x G A*, and t > 0. 

So choosing a^{i) > on A^ (i.e. spins pointing to the z direction at the boundary) 
and also initial condition mo(x) > on A*, Proposition 5.1 shows that m^{t,x) > stays 
bounded away from zero uniformly in time, so is the case for the limiting orbit m(x) on Z^, 
thus D{m) = 0. Our conjecture is again comforted by the following numerical experiments, 
which also show that m(x) depends in a more essential way on the initial conditions than for 
the planar rotator. 
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Fig 10. a: L* 



128, 



1, XY plane 



Fig 10. b: L* 



128, d 



1, YZ plane 



We start with prescribing the spins variables on A'^ as in (4.3), taking a family of loops 
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0, 1, 2 ■ ■ ■ along which cr. 



(cos $ij sin 6i, sin sin 6i, cos 6i), $ij 



27rdj/N, + 



4>Q, and < 9i < 9o, decreasing with i, 6*^;, = on the last loop interacting with A, and 
9q small enough to fit with Belavin-Polyakov conditions. So fixing the precession number 
d = D^{m), we get a "stationary spin wave pattern" on the boundary. Inside A* we choose 
random initial values, |mo(x)| < m/j. 
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Fig 11. a: L* 



128, d 



1, XY plane 



Fig 11. b: L* 



128, d 



1, YZ plane 



We represent here a few sample of (x, y) and (y, z) projections of the field m, which 
yield the following observations. In general, the solution is very sensitive to the choice of 
initial conditions, and many patterns show up, which refiects the moduli in Belavin-Polyakov 
model. For relatively small /3 (e.g. /3 = 5 with L* = 128) spin waves fiuctuate, and can 
take negative values, but the domain where spins point downwards is not sufficiently large 
to start revolving around the sphere. So the degree is D = 0. We still get 2 dimensional 
"vortices" in the (x,y) plane, there are typically 1, 2 or 3 such "vortices" when d = 1, and 
up to 4 when d = 2. Exceptionally, we can also get a 3 dimensional vortex, i.e. xq such that 
m(xo) becomes quite small. Such a m is no longer homotopic to a function on the sphere. 



Increasing /3 generally prevents getting too small values for m, and allows larger negative 
m^. For P = 10 and L* = 128, there are random trials where the winding number D is non 
zero. Thus Fig. 10 is obtained for d = 1 and suggests also D = 1. In Fig. 11 we still have d = 1, 
but D = 0, although 2 large symmetric regions contain negative values of m^. Fig. 12 gives 
an example where D = 1 for d = 2. Nevertheless, we have not observed winding numbers 
D = 2 for = 2. 
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Fig 12. a: L* = 128, (i = 2, XY plane 



Fig 12. b: L* = 128, (i = 2, YZ plane 



Further increasing /3 for a given L* doesn't reveal anything new ; namely, if small temper- 
ature seems to favors long range order and existence of non trivial instantons, it also creates 
stiffness and a need for space. In any case, one should keep in mind that Belavin-Polyakov 
instantons can be reproduced only as /3 ^ oo, and in the thermodynamic limit |A| ^ oo. Of 
course, everything can be again improved through simulated annealing. 
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